Sky localization of complete inspiral-merger-ringdown signals for nonspinning massive 

black hole binaries 
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We investigate the capability of LISA to measure the sky position of equal-mass, nonspinning 
black hole binaries, combining for the first time the entire inspiral-merger-ringdown signal, the 
effect of the LISA orbits, and the complete three-channel LISA response. We consider an ensemble 
of systems near the peak of LISA's sensitivity band, with total rest mass of 2x 10 6 M , aredshift of 
z = 1, and randomly chosen orientations and sky positions. We find median sky localization errors 
of approximately ~ 3 arcminutes. This is comparable to the field of view of powerful electromagnetic 
telescopes, such as the James Webb Space Telescope, that could be used to search for electromagnetic 
signals associated with merging massive black holes. We investigate the way in which parameter 
errors decrease with measurement time, focusing specifically on the additional information provided 
during the merger-ringdown segment of the signal. We find that this information improves all 
parameter estimates directly, rather than through diminishing correlations with any subset of well- 
determined parameters. Although we have employed the baseline LISA design for this study, many 
of our conclusions regarding the information provided by mergers will be applicable to alternative 
mission designs as well. 
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I. INTRODUCTION 

Among the gravitational-wave (GW) sources which the 
proposed Laser Interferometer Space Antenna (LISA) 
would observe, massive black hole binaries (MBHBs) 
stand out in several ways. They have the highest typical 
signal-to-noise ratios (SNRs) and will be observable out 
to redshifts of z ~ 10 or greater JlJ-Q , making MBHBs 
the most distant sources that LISA will observe. Obser- 
vations of MBHBs figure prominently in LISA's science 
goals [3] and may provide particularly rich opportunities 
for multi-messenger observation in conjunction with elec- 
tromagnetic instruments. The potential for such obser- 
vations depends particularly on LISA's ability to localize 
the gravitational-wave sources on the sky and thereby 
facilitate simultaneous or follow-up electromagnetic ob- 
servation Q. 

LISA is an all-sky instrument in the sense that, at any 
one instant, it is sensitive to sources at most points on 
the sky. While this is advantageous from the perspec- 
tive of detection of unknown sources, it is not helpful 
for locating these sources on the sky. To build localiza- 
tion information, LISA takes advantage of the fact that 
MBHBs are observable for months or even years prior to 
merger. During this time, LISA's orbit around the Sun 
introduces both frequency and amplitude modulations to 
the signal that can be used to determine source location 
Q. The majority of studies of LISA's source location 
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abilities [6l-[h5| focus on this mechanism and neglect the 
final merger-ringdown portion of the MBHB waveform, 
which is difficult to model and too brief to experience any 
modulation from LISA's orbit. 

Recent work [l6|, [13] has shown that the merger- 
ringdown portion of the MBHB signal, which is short- 
lived but has more power and occurs at higher- 
frequencies, can achieve sky localization comparable to 
the long inspiral portion without any additional infor- 
mation from LISA's orbital modulation. The mechanism 
for extracting this information is believed to be related 
to the complex response of LISA to gravitational waves 
at the upper end of the measurement band (above ~ 3 
mHz, where finite arm length effects begin to become im- 
portant [l8|). In this paper, we investigate the combined 
effect of both localization mechanisms. Specifically, we 
expand our previous analysis to include the full inspiral 
signal, including orbital effects. We also include a third 
signal channel, T, relevant only at higher frequencies, 
which was excluded from the earlier analysis for techni- 
cal reasons. 

The outline of the paper is as follows: In Sectio n HI1 
we review our methodology, which closely follows [17| . 
Section III Al describes the construction of the inspiral- 
merger-ringdown (IMR) waveform using the "IRS-EOB" 
model Section lll Bl describes our model for the LISA 
instrument, including the response function and noise. 
Finally, Section lll Cl explains how we calculate parameter 
errors using the standard Fisher matrix formalism. 

In Section Mil we examine LISA's ability to localize 
sources, including the improvements that arc gained as 
different levels of realism arc added. Wc find that for an 
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ensemble of systems with total rest mass of 2 x 10 6 M Q 
at a redshift of 2 = 1, observing the IMR signal with 
a stationary detector (and therefore without the benefit 
of orbital modulation) provides better sky localization 
than observing the inspiral signal with an orbiting detec- 
tor. We then consider the observation of a complete IMR 
signal with an orbiting detector and find the final error 
to be only a modest improvement over the IMR signal- 
stationary detector result. Measuring the long-duration 
inspiral is essentially irrelevant to achieving LISA's lo- 
calization potential. Nevertheless, the orbital informa- 
tion contained in the inspiral is still important for local- 
izing signals in advance of merger, which may be crit- 
ical for conducting counterpart searches. We find, in 
agreement with previous results, that LISA can local- 
ize most MBHBs to within the Large Synoptic Survey 
Telescope (LSST) field of view (FOV) several months 
before merger. With the complete IMR signal, ~ 50% 
of MBHBs in our ensemble can be localized into the 
much smaller FOV of the James Webb Space Telescope 
(JWST). This is possible mainly because of the extra 
information provided by the T channel. 

In Section IIV1 wc demonstrate that the merger- 
ringdown (MR) portion of the waveform provides addi- 
tional information beyond the inspiral for all parameters. 
We summarize and conclude in Section [Vj focusing on 
how our results may impact alternative designs for LISA. 



II. METHODOLOGY 

The procedure used to estimate the errors with which 
LISA will measure astrophysical source parameters of 
MBHBs closely follows that of our previous studies 
fl6l [l7| . We will provide a brief overview; those inter- 
ested in more detail should consult [l7| . To assess pa- 
rameter accuracy, we require a model for the emitted 
waveform, a model for the detector response and noise, 
and a method for converting these pieces of information 
into a theoretical limit on the achievable accuracy in mea- 
suring parameters with LISA. 



A. Waveform Model 

We use a complete IMR waveform model tuned to 
match the available numerical simulations for nonspin- 
ning black hole binaries. This model, referred to as the 
IRS-EOB model, uses a conventional effective-one-body 
(EOB) Hamiltonian formalism to model the adiabatic in- 
spiral [20| . The merger-ringdown is modeled using the 
"implicit rotating source" (IRS) formalism, which treats 
the radiation source as a shrinking rigid rotator [l9j |. 

The waveform amplitude is calculated using a flux 
model that is constrained both to be consistent with the 
inspiral flux through 3.5 post-Newtonian (PN) order and 
also to vanish as it approaches the ringdown frequency 
(referred to as "Model 2" and given by Eq. 19 in 



The waveform phase is fit to a physically motivated func- 
tional form (see Eq. 9 in [Hj]). 

Working in geometrized units where G = c = M = 
mi + ni2 = 1, the source model depends only on the 
mass ratio q = mi/rri2 (where mi < 7712) and the spins. 
We generate a waveform time series using the IRS-EOB 
model with a signal cadence corresponding to a quarter 
wavelength at the highest frequency reached by the £ = 4, 
m = +4 harmonics 1 . After the source calculation, wc 
transform the waveform from the source coordinates to 
the solar system barycenter (SSB). The SSB waveform 
hs is given by: 
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2 Ye m are the spherical harmonics of spin- weight 
and hi m is the dimcnsionless strain decomposed 



where 
-2 

in that basis in the source frame. This transformation 
depends on six parameters: the redshifted total system 
mass M = M {1 + z), the luminosity distance Dl, the 
coalescence time t c , and three angles describing the orien- 
tation of the binary. The coalescence time, t Cl is defined 
as the time at which the peak of the waveform reaches 
the solar system barycenter. For the angles, we use the 
inclination t (with the convention that 1 = corresponds 
to the line of sight being coincident with the orbital axis 
of the binary), initial orbital phase <p , and polarization 
phase ip. 

We focus on a specific case that has appeared elsewhere 
in the literature (e. g. @,H|): an equal- mass, nonspinning 
MBHB with rest mass M D = 2 x 10 6 M at z = 1, 
observed for the final year prior to merger. This mass is 
interesting because it lies near the most sensitive range 
of the LISA instrument [1| , merger-tree models predict 
LISA will observe a significant number of events in this 
range [22j , and excellent observational evidence exists for 
black holes of this size [23Tj25| . 

The sky location of the binary, described by the ecliptic 
latitude /3 and longitude A, is applied by the instrument 
response, which we discuss below. We use the vector 
A a = (InM, \uDl, P, A, l, 4> , ip, t c ) to denote the com- 
plete set of variable parameters. We exclude the mass 
ratio q from A a , as it is not varied when computing pa- 
rameter uncertainties, a procedure consistent with that 
used in [lE [H HI . 



B. Instrument Model 

The instrument model includes the transformation of 
the SSB waveform hs to the LISA instrument outputs as 
well as a model of the instrument noise. The instrument 



1 We oversample to accommodate interpolation when applying the 
response. 
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outputs are the {A, E, T} TDI combinations described in 
17 and based on the orthogonal {A, E, T} combination 
27 . The response is computed by employing Synthetic 
LISA [28| to generate the {X, Y, Z} TDI combinations 
[29| and then converting to {^4, E, T} using 
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We use an analytic model of the noise in the {^4, E ) T} 
TDI channels, obtained following the procedure in j30| . 
The expressions for the one-sided spectral densities are 



S A<B = 2sin 2 ($) [2(3 + 2 cos($) 
S f = 8sin 2 ($)sin 2 ($/2) [4sin 2 ($/2)S' 1 
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where $ = 2irfL/c, f is the frequency, L is the arm 
length, and S pm and S op are the one-sided spectral den- 
sities of the proof-mass acceleration and optical path- 
length noises, 
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In our prior studies [lfj, [17J , we neglected the T chan- 
nel due to inconsistencies that were observed between the 
analytic noise model for Sf in ([3]) and the noise spectrum 
obtained from a time-domain simulation of Sf performed 
with Synthetic LISA. We have now verified that the 
two models are consistent for / > 0.5 mHz. Where noted, 
we include contributions from the T channel at these fre- 
quencies. We expect this to be a reasonably accurate ap- 
proach since the response of the T channel should vanish 
at low frequencies. 

For the foreground of gravitational waves from unre- 
solved compact binaries, we use the model developed in 
[3lj ] . by adding 



5 ga i = [4$sin($)] 2 5 conf 



(5) 



to the expressions for e ■ ^conf is a signal confusion 
noise estimate for a related measurement variable, taken 
from Eq. 14 of [HJ. No gravitational- wave foreground is 
added to the T channel since S^onf "C Sf at the frequen- 
cies where T is included. 



C. Error estimation 

To approximate the achievable measurement accuracy 
for LISA, we apply the Fisher matrix formalism. While 



not as robust or informative as more sophisticated meth- 
ods such as Markov Chain Monte Carlo [32j , the compu- 
tational efficiency of the Fisher matrix approach is well 
suited for this type of performance estimation work. Po- 
tential pitfalls have been discussed in the literature (33|. 
The Fisher information matrix is defined as 
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where ft. is a channel of TDI strain, a and b are parameter 
indices, and (. . . | . . .) denotes a frequency-domain inner 
product weighted by the inverse noise spectral density 
S n (f) of the TDI output n = A, E, f: 
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Here /i ow is the low-frequency cutoff and /high is the half 
the sampling frequency described in Section III Al The 
individual Fisher matrices for each TDI channel to be 
included in the analysis are summed to form the total 
Fisher matrix. Throughout this work, we calculate the 
parameter derivatives in ([6]) using second-order centered 
differencing (rather than one-sided as in [13]), with a frac- 
tional step size e a = AA a . We set A = 1CT 4 for the 
coalescence time and A = 10 ~ 6 for all other parameters. 

The covariance matrix, E ah , is well approximated by 
the inverse of the Fisher matrix for signals with large 
SNR: 



(8) 



so that a a = VS aa is the standard deviation of parameter 
A a . Because we use lnM and lnD^ as parameters, the 
resulting uncertainties are fractional: 



lnM A/ / ji /T 

a = a /M . 
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The uncertainty of the other dimensionful parameter, t c , 
is not expressed fractionally, as one is generally interested 
in the absolute timing error, but we emphasize that one is 
free to make either choice. We note that a poor choice of 
units for t c could cause the Fisher matrix to have a large 
dynamic range and be ill conditioned for inversion. We 
find that expressing t c in seconds yields computationally 
invcrtible Fisher matrices for the cases we consider. 



III. SKY LOCALIZATION 

As mentioned in Section U the localization of a par- 
ticular MBHB source on the sky is critical for enabling 
coordinated observations with electromagnetic observa- 
tories. One goal might be to locate the host galaxy in 
which the MBHB merger occurred. If the galaxy can 
be identified and the redshift obtained, it can be com- 
bined with the luminosity distance determined from the 
gravitational waveform for use as a cosmological probe 
0, HI| • The simplest way to identify a host galaxy would 
be for LISA to generate an error volume that only con- 
tained a single galaxy. Unfortunately, even the most op- 
timistic estimates of LISA's localization ability suggest 
that this will not generally be possible @. Identifying 
the host galaxy within LISA's error volume will then re- 
quire other information, perhaps the detection of a tran- 
sient electromagnetic signal associated with the merger 
itself. The nature of this possible signal is a subject of ac- 
tive research [35l [36j . and it is unclear if it exists, when 
it might peak relative to the GW signal, and what its 
spectral content will be. 

The observational effort needed to search for an EM 
counterpart to an MBHB merger will depend on the ra- 
tio between LISA's sky localization error box and the 
field of view (FOV) of the EM asset used to perform the 
search. Roughly speaking, this ratio will determine the 
number of telescope pointings that are needed to cover 
the LISA error box. While some counterpart scenarios 
involve afterglows that persist long after coalescence, oth- 
ers involve signals that peak at or before coalescence and 
may fade quickly. For these signals, it is important to 
have some advance warning as to where on the sky (and 
when) the merger will occur, allowing for appropriate 
EM resources to be scheduled and pointed in advance. 
The sky localization information produced by LISA's or- 
bit accumulates with time over the months or years that 
the MBHB is observed, allowing for such advance notices 
to be issued. 

Figure [T] summarizes various estimates for the devel- 
opment of precision in ecliptic latitude over time for a 
LISA observation of two equal-mass, nonspinning black 
holes with a combined rest mass M Q = 2 x 10 6 M Q at 
z = 1. All of the curves show the median results from 
an ensemble of 100 systems with random sky locations 
and orientations. Except where noted, the low-frequency 
limit of the inner product ([7]), /i OW! was set to 0.03 mHz. 
The precise placement of this cutoff and the model of the 



instrument noise between the cutoff and the lower limit 
of the "official" LISA measurement band at 0.1 mHz can 
have a significant effect on the results [!, [|[ . The curves 
are all shown as functions of time before i c , defined as 
the peak of the waveform. In addition, we also show to 
the right the final error achieved with the complete MR 
signal. 

The purple dot-dashed curve with square markers uti- 
lizes PN waveforms and a simplified model of the LISA 
instrument response known as the low-frequency approx- 
imation (LFA). This approximation produces two chan- 
nels of information which can be compared to the A and 
E TDI channels discussed in Section Hi B I The estimate 
is generated by the same code used in [ill. [l3| , with the 
exception that here we ignore the effects of spin preces- 
sion. We include it primarily for comparison with our 
own method described in Section [TTJ Due to the break- 
down of the post-Newtonian approximation, this curve is 
truncated at the time corresponding to the Schwarzschild 
innermost stable circular orbit (ISCO). Therefore, this 
estimate is sensitive to sky localization information en- 
coded in the waveform by LISA's orbital motion but is 
not sensitive to localization information resulting from 
inclusion of the MR waveform or from the high-frequency 
response of the instrument. 

The red dotted curve with diamond markers, based 
on our previous work [l7j . uses the IRS-EOB waveform 
model, as described in Section Hi Al and the Synthetic 
LISA instrument response described in Section al Bl Due 
to the aforementioned problem (now solved) of spurious 
information in the T channel, only the A and E TDI 
channels are included. In addition, the instrument re- 
sponse uses a stationary, non-orbiting constellation lying 
in the ecliptic plane instead of the usual LISA orbits. 
This estimate is therefore only sensitive to the sky lo- 
calization information provided by the interaction of the 
MR signal with the high-frequency response of the de- 
tector. Lacking the orbital contribution, we can safely 
remove a large portion of the inspiral by setting the low- 
frequency limit to /i ow = 0.1 mHz for this single estimate 
only. This helps reduce computation time. 

For times before the ISCO, the errors in the orbits-only 
estimate (purple, dot-dashed, squares) are far less than 
those in the stationary-LISA estimate (red, dotted, dia- 
monds), suggesting that the orbital mechanism for sky lo- 
calization dominates over these time periods. The dotted 
curve only "catches up" to the dot-dashed curve ^100 
seconds before the merger. However, by the time the 
full signal is considered, the stationary-LISA errors have 
become ~ 4 times smaller than the orbits-only errors at 
ISCO. In other words, for nonspinning, equal-mass bina- 
ries in this mass range, observing the MR signal alone 
gives better localization information than observing the 
entire inspiral. 

The remaining curves are new results, estimates which 
include both the complete IMR (IRS-EOB) waveform 
and a realistic, orbiting (Synthetic LISA) model of 
the instrument. In addition, the low frequency cutoff has 
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FIG. 1: Comparison of latitude uncertainty estimates for an equal-mass, nonspinningMBHB with total mass 2 x 10 6 Mq 
and redshift 2=1. The purple dot-dashed curve with square markers (based on [jj], [ll|) uses a post-Newtonian inspiral (I) 
waveform and a low-frequency approximation (LFA) to the detector response that includes LISA's orbit. The red dotted curve 
with diamond markers (based on [13]) uses the complete inspiral-merger-ringdown (IMR) waveform and Synthetic LISA (SL) 
for the response function but neglects LISA's orbit, all frequencies below 0.1 mHz, and the T channel. The green dashed curve 
with triangle markers uses the IMR waveform, the SL response, and orbital effects, but neglects the T channel. The blue solid 
curve with circle markers includes all effects: the IMR waveform, the SL response with orbits, and all three TDI channels. The 
FOVs of LSST and JWST, as well as the location of the Schwarzschild ISCO, are shown for reference. 



been restored to 0.03 mHz in order to take advantage of 
the orbital modulation. The green dashed curve with tri- 
angle markers includes only the A and E TDI channels, 
while the blue solid curve with circle markers includes 
the T channel as well. 

The A , E estimate (green, dashed, triangles) tracks 
the orbits-only estimate (purple, dot-dashed, squares) at 
times well before merger, thus confirming that our calcu- 
lation is consistent with that of Lang and Hughes [Tl[ [l3| 
in the inspiral regime. It also tracks the stationary- 
LISA estimate (red, dotted, diamonds) as the MBHB 
approaches merger. When the entire IMR waveform is 
included, the orbiting and non-orbiting A , E estimates 
differ by only ~ 10%. The dashed curve contains both 
mechanisms for localizing a source (orbits and high- 
frequency detector response to merger- ringdown), but 
it ultimately does not do significantly better than the 
stationary-LISA response alone. Our ability to localize 
sources well in advance of merger depends critically on 
the orbital information. However, for the final achiev- 
able sky location accuracy in these systems, the inspiral 
is essentially irrelevant and only the MR signal need be 
taken into account. 

The estimate with all three TDI channels (blue, solid, 
circles) begins to diverge from the orbits-only estimate 
(purple, dot-dashed, squares) at earlier times than the 
A , E estimate and reaches a value at coalescence that is 
about half that of the A , E estimate. This represents the 
additional contribution of the T channel, which is con- 
centrated at high frequencies. It is neglected implicitly in 
the orbits-only estimate through the use of the LFA and 



explicitly in the other two estimates. While all three TDI 
channels are available in a fully operational, 6-link LISA, 
only two channels are available if one link is down j38j . 
Using the green dashed line as a proxy for a 5-link LISA 
and comparing to the blue solid line, we can see roughly 
how much the sky position determination is affected by 
a non-optimal detector. 

Figure [T] also includes the FOV of two representative 
EM instruments that could be used for locating coun- 
terparts. The Large Synoptic Survey Telescope (LSST) 
[39l | has a FOV of 9.6 deg 2 , which we have modeled as a 
square aperture ~ 186 arcmin on a side, and is designed 
to rapidly and repeatedly scan large areas of the sky. 
The figure demonstrates that LISA will localize a typi- 
cal equal-mass, nonspinning MBHB source of this mass 
and redshift within a single LSST field several months in 
advance of merger. This will allow multiple exposures of 
the target area to be made throughout the event in an 
effort to detect an optical transient. 

If a transient signal is not bright enough to be visible 
with a survey telescope such as LSST, then a more sensi- 
tive instrument with a narrower field of view may be re- 
quired. In Fig.[T]we plot the 2.2 arcmin x 4.4 arcmin FOV 
of the NIRCam instrument on the James Webb Space 
Telescope (JWST) [4(| as an example of such an instru- 
ment. Our most complete estimate of LISA's localization 
capabilities (IMR waveform, orbiting detector with Syn- 
thetic LISA response, all TDI channels) has a median 
latitude error of 2.6 arcmin. This suggests that a single 
JWST exposure might be sufficient to cover the LISA 
error box for this system. In reality, a few exposures 
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FIG. 2: Latitude uncertainty, cr^, for a 2 x 10 6 M Q system at z = 1. The center line (with individual time samples marked) 
denotes the median error, while the shaded region depicts the middle three quintiles (20 th — 80" 1 percentile) of errors. Both 
plots correspond to IMR signals with an orbiting LISA constellation. However, panel (a) shows results for just the A and E 
TDI channels, while panel (b) also includes the T channel. 




FIG. 3: Same as Fig. [2] but for uncertainty in longitude, cr"* 



may be required to account for the shape of the LISA er- 
ror box, its orientation with respect to the JWST FOV, 
and any operational constraints of JWST. The results in 
Fig. Q] suggest that unlike the situation with LSST, in 
which a number of pre-merger images could be obtained, 
a coordinated observation with JWST (or any other in- 
strument with a similar FOV) will have to be a follow-up, 
as LISA requires the complete IMR waveform to provide 
precise localization information. For example, at ISCO, 
the best estimate for median latitude error exceeds the 
JWST FOV by roughly half an order of magnitude. 

All of the curves in Fig. [T] represent the median LISA 
latitude error for an ensemble of MBHB mergers. Fig- 
ures [5] and [3] also show the distribution of both lati- 
tude and longitude errors, based on estimates with the 
complete IMR waveform and LISA's orbital modulation. 
Separate panels show results including just A and E and 
results including all three TDI channels. In each panel, 
the solid line represents the median values (equivalent to 
what is plotted in Fig. [T]), while the shaded area shows 
the region spanned by the middle three quintiles of the 
distribution. The left panels (A, E only) show that the 
final LISA error box will lie within the FOV of JWST for 
systems in the 80 th percentile of sky localization preci- 
sion. That is, both latitude and longitude are measured 
to within ~ 4.4 arcmin or better for the top 20% of sys- 
tems. When T is included, the situation is better still, 
with a full 50% of binaries in our ensemble fitting within 



the JWST FOV 2 . These figures show just how important 
T is for computing final localization errors. Not only do 
the medians drop for both parameters, but the width of 
the distribution decreases as well. Although a 5-link and 
6-link LISA are indistinguishable at early times, com- 
plete functionality can have important effects for the full 
signal. 

For completeness, in Fig. 0] we show errors in lumi- 
nosity distance, cj Dl /Dl, for IMR signals with orbital 
modulation. Distance is the third axis of the 3D local- 
ization "voxel." Given an assumed cosmological model, a 
measurement of distance is equivalent to a measurement 
of redshift. This approximate rcdshift can help decide be- 
tween potential host galaxies located within the 2D error 
box. We see that the median distance error drops below 
a percent within the last day, reaching a minimum of less 
than 0.1%. Unfortunately, this "intrinsic LISA error" 
is not the limiting uncertainty for measuring luminosity 
distance. Weak lensing from matter between the source 
and the detector will slightly magnify or demagnify the 
GW signal in an unknown way, adding significant error 



2 These assertions are somewhat crude: The plot does not actually 
show whether individual binaries meet the criterion. It is possi- 
ble that some systems are well localized in latitude but poorly 
localized in longitude, and vice-versa. A more careful check of 
the raw data verifies the numbers in the text. 
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FIG. 4: Same as Fig. [2] but for fractional luminosity distance uncertainty, a° L / Dl . An estimate of the weak lensing limit 
(WL) [37J at z = 1 is shown for reference. Figs. [2^4] collectively show the evolution of the 3D error voxel. 



to the distance measurement. The line marked "WL" on 
the figure is an estimate of <t Dl /Dl ~ 4.4% from 
weak lensing at z = 1. The intrinsic LISA error crosses 
the WL error over a month before merger; shortly af- 
ter that time, the information on Dl provided by LISA 
becomes irrelevant unless weak lensing error can be mit- 
igated significantly |4lj . 



IV. CONTRIBUTION OF MERGERS TO 
PARAMETER ESTIMATION 



The results presented in the preceding section clearly 
demonstrate that observing the MR portion of an MBHB 
signal can provide a significant amount of sky localization 
information, despite the fact that there is little or no 
modulation imposed by the LISA orbits. In this section, 
we explore the nature of this information. Throughout 
this section, we always use the complete IMR waveform 
with orbital modulation and all three TDI channels. 

One question of interest is whether the information 
provided by the MR signal is tightly coupled to just one 
or two of the elements in A a or if it generally improves 
the knowledge of all parameters. Figure [5] shows the time 
evolution of the median eigenvalues of the Fisher matrix 
for the ensemble of systems studied in Section IIIII At 
each time step, r a f, was calculated using the A, E, and 
T observables for every case in the ensemble using the 
procedure described in Section [TTJ These matrices were 
used to compute the eigenvalues for every case in the 
ensemble, and the median values are plotted. A similar 
procedure was used to compute the median eigenvectors 
corresponding to these eigenvalues. 

As © indicates, the Fisher matrix is a measure of 
inverse error or "information." The increasing trends in 
the curves in Fig. [5] represent information accumulating 
in the LISA observables as the signal evolves, just as 
the decreasing trends in Figs. Q3H] represent a decrease 
in measurement error. The largest eigenvalue in Fig. [5] 
(solid curve) corresponds to an eigenvector that mostly 
consists of one parameter, InM. The majority of the 
information in this eigenvalue accumulates early in the 
inspiral phase of the signal. From ISCO through the end 
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FIG. 5: Evolution of the median Fisher matrix eigenvalues 
with time for the inspiral, merger, and ringdown of an ensem- 
ble of equal-mass, nonspinning MBHBs with M = 2 x 10 6 Mq 
at z = 1 observed by an orbiting, six-link LISA via the A, E, 
and T TDI channels. 



of the signal, the eigenvalue roughly doubles, with most of 
the increase occurring very near the time of coalescence. 
This is consistent with previous results that show that 
InM is primarily determined during the inspiral phase 
[4^ and that the MR signal provides little additional 
information (l7j . By comparison, the increase in SNR 2 , 
which one might naively expect to be proportional to 
Fisher matrix eigenvalues, is ~ 100 over this period. 

The band of eigenvalues in the middle of Fig. [5] (dashed 
curves) represent a mix of the five angular parameters in 
A a and lnD^. The increase in these eigenvalues from 
ISCO through ringdown is much more significant than 
that for InM, ranging from two to three orders of mag- 
nitude. 

The smallest eigenvalue (dotted curve) corresponds to 
an eigenvector that mostly consists of t c . The small size 
of this eigenvalue is partly due to the intrinsic scale of the 
parameter, which was the only one in A a to be expressed 
in physical units (seconds). It increases by a factor of 
~ 30 between ISCO and the end of ringdown. 

The fact that all of the eigenvalues increase post-ISCO 



8 



indicates that the MR portion of the waveform, coupled 
with the high-frequency response of the LISA instrument, 
provides information in all directions of parameter space. 
We might have expected that the merger, with its pro- 
nounced peak, would provide a very accurate measure- 
ment of t c and only affect other parameters by decreas- 
ing their degeneracy with t c . However, Figure [S] demon- 
strates that this is not the case. If we were to assume 
that sky localization results primarily from the modu- 
lation imposed by LISA's orbit, then we would expect 
that the merger cannot provide any intrinsic information 
about sky position, since the LISA constellation remains 
effectively stationary during merger. Figure [5] demon- 
strates that the interaction between the MR signal and 
the complex response of the LISA instrument provides 
information about every parameter independently. This 
is in addition to any contribution from degeneracy break- 
ing which may occur. 

V. SUMMARY AND CONCLUSIONS 

We have presented the first estimates of LISA's ability 
to measure the source parameters of nonspinning, equal- 
mass MBHBs that include a complete inspiral-merger- 
ringdown waveform, realistic LISA orbits, and the full 
high-frequency response of the LISA instrument. For sys- 
tems with M a = 2 x 10 6 M Q at z = 1, we find that LISA 
can locate the source on the sky with an error of ~ 3 
arcminutes in the median case at the end of the MBHB 
event. We have also studied how the sky localization in- 
formation for these systems evolves with time and find 
that, at ISCO, the sky localization errors are typically 
tens of arcminutes. 

Inspired by the potential for multi-messenger astron- 
omy enabled by combined observations of MBHB mergers 
with gravitational-wave detectors and electromagnetic 
telescopes, we have compared the source-localization ca- 
pabilities of LISA with the FOVs of two future tele- 
scopes: LSST and JWST. We find that most systems 
in our ensemble are located within the LSST FOV sev- 
eral months prior to coalescence. Localization within the 
smaller FOV of JWST requires information from the MR 
waveform. 

While the results in this paper focus on a particular 
instrument (LISA) and a particular source class (equal- 
mass, nonspinning, M a = 2 x 10 6 M , z = 1) they can 
provide insight for other scenarios as well. For example, 
proposals for lower-cost variants of LISA often involve a 
constellation with shorter arm lengths [43| . This has the 
general effect of shifting the response of the instrument to 
higher frequencies. We have asserted that the improve- 



ment in parameter estimation that we observe post-ISCO 
is due in part to the complexity of the high-frequency re- 
sponse of the instrument. If this is true, shortening the 
arms could have a significant impact on the ability to 
localize systems in the mass range studied. Specifically, 
the improvement in sky localization from observing the 
merger-ringdown of a system with M a ~ 10 6 M Q may 
not be as dramatic as that presented in Section Hill At 
the same time, the decrease in strain sensitivity associ- 
ated with shorter arms would reduce the duration of the 
inspiral signal that is observed, further degrading mea- 
surement performance for these systems. For systems 
with larger M a , the situation would be worse. 

For lower-mass systems, shortening the arms (if the in- 
terfcrometric measurement system is not simultaneously 
degraded) may actually increase the benefits of observ- 
ing the merger by increasing the SNR for the merger- 
ringdown portion of the signal. Since these lighter sys- 
tems have higher merger frequencies, they will still ex- 
perience a complex response in the detector, even with 
shorter arms. Determining whether this effect can over- 
come the lower intrinsic SNR of less massive systems will 
require a more detailed calculation. 

Finally, we note that while the addition of the merger- 
ringdown signal and a complete detector response is an 
important step towards realizing the full potential of 
LISA and other space-based gravitational-wave detec- 
tors, we still do not have a complete picture. Other po- 
tential sources of information in the MBHB waveform, 
particularly spin and the associated precession effects, 
have been demonstrated to significantly improve LISA's 
localization ability in the inspiral-only case Q. It re- 
mains to be seen what impact combining spin precession 
and merger will have on LISA sky localization estimates. 
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